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The recently completed High Altitude Water Cherenkov (HAWC) gamma-ray observatory has 
been taking data with a partial array for more than one year and is now operating with >95% 
duty cycle in its full configuration. With an instantaneous field of view of 2 sr, two-thirds of 
the sky is surveyed every day at gamma-ray energies between approximately 100 GeV and 100 
TeV. Any source location in the field of view can be monitored each day, with an exposure of 
up to ~ 6 hours. These unprecedented observational capabilities allow us to continuously scan 
the highly variable extra-galactic gamma-ray sky. By monitoring the flaring behavior of Active 
Galactic Nuclei we aim to significantly increase the observational data base for characterizing 
particle acceleration mechanisms in these sources and for studying cosmological properties like 
the extra-galactic background light. In this work we present first studies of data taken between 
June 2013 and July 2014 with a partial array configuration. Flux light curves, binned in week- 
long intervals, for the TeV-emitting blazars Markarian 421 and 501 are discussed with respect to 
indications of flaring states and we highlight coincident multi-wavelength observations. Results 
for both sources show indications of gamma-ray flare observations and demonstrate that a water 
Cherenkov detector can monitor TeV-scale variability of extra-galactic sources on weekly time 
scales. The analysis methods presented here can provide daily flux measurements with a mini¬ 
mum time interval of one transit and will be applied to new data from the completed HAWC array 
for monitoring of blazars and other transients. 
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1. Introduction 

Active Galactic Nuclei (AGN) are the most populated source class of very high energy (VHE; 
> 100 GeV) gamma-ray emission. 1 Those AGN with jets oriented close to the line of sight 
are called blazars and among them are some of the brightest VHE sources. VHE fluxes from 
blazars are often highly variable, and periods of increased activity, called flares from here on, have 
been observed with fluxes increasing by an order of magnitude [|I]] and with variability on time 
scales from months down to minutes [Q, j31. It is poorly understood where along the jet these 
gamma rays are emitted and if the underlying processes that accelerate particles to TeV or higher 
energies are dominated by populations of leptonic or hadronic particles. Only a more detailed 
coverage of flares over a wide range of wavelengths can constrain such models as in [14!]. Most of 
the observations at VHE energies have been performed with Imaging Air Cherenkov Telescopes 
(IACTs) that can only operate during clear, dark nights and normally monitor only one source in 
the field of view at any time. Their capacity for long term studies of the frequency of blazar flaring 
states is therefore limited. Furthermore, extended monitoring that is not biased by following-up 
on flaring alerts at other wavelengths can help to study the possibility of flares limited to TeV 
energies, so-called orphan flares, which might help to identify hadronic acceleration processes [5]. 
Studies of the flux variability in AGN and their underlying acceleration mechanisms will greatly 
profit from unbiased monitoring observations with high duty cycles that are made possible with the 
HAWC Observatory. In addition, bright flares offer a unique opportunity to constrain or measure 
the extra-galactic background light [6] and inter-galactic magnetic fields [|7j], 

2. Analysis method 

2.1 The HAWC Observatory 

The High Altitude Water Cherenkov (HAWC) Observatory is located at an altitude of 4,100 m 
above sea level on the slope of the Sierra Negra volcano (19.0°N 97.3°W) in the state of Puebla, 
Mexico. The detector is optimized for measuring extensive air showers induced by gamma rays 
with energies between approximately 100 GeV and 100 TeV. Completed in March 2015, HAWC 
is an array of 300 Water Cherenkov Detectors (WCDs). Each WCD is composed of a large steel 
tank holding ~ 200,000 liters of purified water in a light-proof bladder. Four photo-multiplier 
tubes (PMTs), located at the bottom, can detect the Cherenkov light from the particles in an air 
shower passing through the array. By analyzing the recorded light intensity distribution as a proxy 
for the particle density and fitting the plane of the shower front from the photon arrival times it 
is possible to reconstruct the size and incident direction of an air shower event. While most of 
the air showers recorded are produced by hadronic primaries, this background can be significantly 
reduced by identifying large charge deposits outside the core regions indicative of muons, which 
are not expected in gamma-ray showers. More details about the reconstruction, data taking and 
hadron rejection are discussed elsewhere in these proceedings |[sj jh ]. 

The modular design of HAWC made it possible to start continuous scientific operation before the 
completion of the whole detector. The results presented here are from data taken between June 13, 

1 See list at http://tevcat.uchicago.edu . 
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2013, when HAWC was operating with 307 active PMTs in 90 WCDs, and July 9, 2014, by which 
time the array had grown to include 495 active PMTs in 135 WCDs. During the whole period 
the detector was operating with a duty cycle ~ 90 % with down time mostly due to infrequent 
maintenance activities that required a complete shutdown. 


2.2 Data reduction 


After calibration and reconstmction, the air shower events are sorted into 10 analysis bin, 
primarily defined by the fraction of PMTs that measured a signal in the event. The lower threshold 
of the first analysis bin is 4.5% of all active channels, thus requiring ~ 20 PMTs with photon hits 
in the same event. As a measurement of shower size, the number of PMTs with signals serves 
both as a rough proxy of the energy of the primary particle and as an event quality parameter, 
since larger events with more signals will on average have a more precise angular reconstruction. 
For each analysis bin, additional cuts on several event parameters reduce the number of hadronic 


background events, as mentioned in section [24 


Due to the increasing absorption of secondary particles in the atmosphere, the sensitivity to gamma 
rays is a function of zenith angle. For the time-dependent analysis in this contribution we use 
the duration of one source transit over HAWC as the smallest time interval under consideration, 
thus simplifying this dependence to be a function of the source declination. A HAWC analysis 
targeted for blazars flares with durations of less than one transit is described elsewhere in these 


proceedings [10]. Since the gamma-ray efficiency of HAWC for zenith angles larger than 45° 
is negligible, a source transit can be considered to last at most ~ 6 hours. In order to store the 
data in time intervals that cover only one transit, the reconstructed events are sorted into sidereal 
days, starting at midnight local sidereal time at the HAWC longitude of 97.3°W. This choice leads 
to transits being split in two for sources with right ascension (RA) < 3 h or > 21 h, which are 
analyzed with a separate binning of the data starting at noon sidereal time. 

For each sidereal day and each of the ten analysis bins, a sky map of event counts is produced by 
populating pixels with an average spacing of ~ 0.1°. In general, these maps are still dominated by 
hadronic background events, and direct integration JIH], averaging counts over two hours around 
any location, yields an estimate of this background. The estimated background counts in each 
pixel are stored in a second map with the same grid structure, so that a difference of the two maps 
provides event excess counts. 


2.3 Likelihood analysis 

To analyze the combination of the 10 analysis bins for a particular region of the sky we employ 
a likelihood ratio test, a standard technique in particle physics [[12|]. The software implementation 
for performing this test on a given set of HAWC data maps interfaces directly with the general 
HAWC software framework and is called Likelihood Fitting Framework (LiFF), described in de¬ 


tail elsewhere in these proceedings Q13p. In this contribution, the only signal model considered 
represents a gama-ray point source with the differential flux energy spectrum described by a power 
law with index s and an optional exponential cut-off c: 
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In LiFF, this input flux is convoluted with a detector response function that includes the point spread 
function and efficiency of triggers and cuts, depending on primary energy and incident angle. For 
one source transit over HAWC, the signal hypothesis contributions as a function of zenith angle are 
summed and yield the expected number of events Sb, p per analysis bin b and pixel p in the map for 
one sidereal day. In cases were the coverage of a source transits is interrupted, for example due to 
a detector shutdown, the lost signal fraction compared to a full transit is calculated by excluding 
the gap period from the integration over zenith angle and this fraction is then used to reduce the 
expected event count accordingly. The likelihood of a hypothesis defined by {Fb,s,c} given an 
observation N of numbers of events in all bins and pixels is expressed via the likelihood function 

C-s ({fodT }is)=nnm. ,p,b. p ) , (2-2) 

V bp 


where P(Nb,p,kb :P ) is the Poisson distribution for a mean expectation A b, p = Sb, P +#/>,/» he. the 
sum of the expected signal and the number of background events estimated from data for analysis 
bin b and pixel p. In the likelihood ratio test, the result of eq. ( |2.2| ) is divided by the likelihood 
value Cb for a background-only assumption (S/, :/ , = 0). In the software implementation, this ratio is 
expressed as a difference A In £ of the logarithms of the two likelihood values and the test statistic 
is defined as TS = 2Aln£ . TS is then numerically maximized by iteratively changing the input 
parameters, yielding those values that have the highest likelihood of describing the observed data 
for the point source model assumption. 

In the analysis presented here, the spectral index s and cut-off value c in eq. ( fP| ) were kept constant 
for each target source and only the normalization Fo was left free in the likelihood fit. This is rea¬ 
sonable since even the one-year time-integrated data has a limited sensitivity to spectral parameters 
due to small statistics. A first look at AGN spectra with HAWC is discussed in JT4|]. 


2.4 Monitoring strategy 


Given the limited sensitivity of the partial HAWC detector, the results of a likelihood fit to 
the data from one transit exhibits large uncertainties in this early analysis. The Crab Nebula, the 
brightest steady source seen by HAWC, yields an average significance of less than 1 a per transit. 
For this reason, the monitoring presented in this contribution is performed for time intervals of 
7 transits, yielding a total of n t = 56 light curve bins. It is therefore insensitive to flux variability 
on shorter time scales, but provides more stable result for a first analysis of the blazar monitoring 
potential of HAWC. A light curve for the Crab Nebula, derived with the methods described here, 


is shown elsewhere in these proceedings | |l 5| [ and is consistent with a constant flux over time. 

The integrated life time of the data in the period June 13, 2013, to July 9, 2014, is equivalent to 
approximately 350 transits, varying slightly (~ 5 days) for different source right ascensions due to 
small variations in the exposure. Among the locations that surpass a pre-trial significance in the 
time-integrated analysis, only two are identified as known extra-galactic sources, the two VHE- 
emitting blazars Markarian 421 and Markarian 501. The monitoring results for these two sources 
are discussed in the following section. 

The Crab Units on the right axes in the light curve figures in section || are defined by dividing the 
measured flux values by the average flux of the Crab nebula as observed in HAWC, derived by 
assuming a simple power law spectrum with s = 2.63 [p~6|] and the flux normalization result for the 
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Days after 2013/06/13 (MJD 56456.0) 


Figure 1: Flux light curve of Markarian 421 for the time between June 13, 2013, and July 9, 2014, in 
intervals of 7 days, with horizontal error bars indicating the time span between the start of the first and the 
end of the last transit over HAWC. Integrated fluxes above 1 TeV for a differential spectrum with power law 
index 2.2 and exponential cut-off at 5 TeV are shown in photons per cm 2 s 2 on the left axis and divided by 
the average Crab flux observed by HAWC on the right axis. The statistical 1 and 2 a flux uncertainties of 
the individual measurements are displayed as orange/yellow bars around the weighted average flux (dashed 
line). The large uncertainties in the flux around day ~ 90 are due to a period of construction and maintenance 
in September 2013 during which HAWC was shut down during day time for several days in a row, creating 
large gaps in the transit coverage. The blue band indicates the period over which a strong increase in X-ray 
flux was observed by MAXI [ 221. 


data for the whole period. In general, we estimate the systematic uncertainty of the HAWC fluxes 
to be ~ 40% which is expected to be reduced with further progress on the data analysis [|]~5|]. 

3. Results 

3.1 Markarian 421 

Markarian 421 (Mrk 421) is a BL Lacertae type blazar with a redshift of z ~ 0.03 that has 
been well studied by many VHE gamma-ray observatories and is known to exhibit a high degree of 
variability [Jl7j]. For this analysis, we use a simple power law signal hypothesis with index 5 = 2.2 
and exponential cut-off at c = 5 TeV, based on typical established values, for example [JT^, |l9j ]. 
The results from the likelihood fit of flux normalizations in each time interval are displayed as 
integrated fluxes in the light curve shown in Fig. [[} 

For the period considered here, the data base of The Astronomer’s Telegram 2 lists only one report 
of increased activity in Mrk 421, an X-ray observation at 2-4 keV by the MAXI/GSC instrument 
during August 2013 [p2fl. The corresponding blue band in Fig. [I] roughly visualizes the length of 
that period of X-ray brightening. The highest flux value obtained in the HAWC analysis coincides 
with the onset of the MAXI observation, several days before their reported peak X-ray flux at 
Modified Julian Date (MJD) 56524. The chance probability of the HAWC measurment alone to 

2 http://www.astronomerstelegram.org 
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Figure 2: Flux light curve of Markarian 501 for the time between June 13, 2013, and July 9, 2014, in 
intervals of 7 days, with horizontal error bars indicating the time span between the start of the first and the 
end of the last transit over HAWC. Integrated fluxes above 1 TeV for a differential spectrum with power law 
index 2.8 are shown in photons per cm 2 s 2 on the left axis and divided by the average Crab flux observed 
by HAWC on the right axis. The statistical 1 and 2 a flux uncertainties of the individual measurements are 
displayed as orange/yellow bars around the weighted average flux (dashed line). The large uncertainties in 
the flux around day ~ 90 are due to a period of construction and maintenance in September 2013 during 
which HAWC was shut down during day time for several days in a row, creating large gaps in the transit 
coverage. The blue/green bands indicating increased X-ray/gamma-ray fluxes observed by other instruments 
are based on publicly available data and discussed and referenced in the text. 


occur as a fluctuation, based on the value of the excess compared to a constant average flux and 
a trial factor for 56 time bins, is 1.4 • 10 5 (4.3a). This is a strong indication that both HAWC and 
MAXI data show the same flaring event, though we have not yet developed a statistical method to 
determine the significance of this multi-wavelength correlation. 


3.2 Markarian 501 


Markarian 501 (Mrk 501) is another BL Lacertae type blazar with a redshift of z ~ 0.03 and 
with a history of known VHE gamma-ray flares, see for example [JT7j] , For the HAWC analysis, 
we used a fixed spectral shape with index s = 2.8 and no cut-off (c = °o) in eq. [24], corresponding 
to typical values [|20j]. The HAWC flux light curve results are shown in Fig. f2j. The highest flux 
value in the HAWC analysis is observed during a week in July 2013 that largely overlaps with a 
period of increased flux observed by the Swift/BAT hard X-ray transient monitor in the 15-50 keV 


band, reported in [|23|]. The left vertical blue band visualizes the approximate extent of the X-ray 
observation that revealed a peak flux at MJD 56495, in the right half of the band. Our a-posteriori 
analysis cannot quantify the probability for this coincidence but seems to indicate that the HAWC 
observation shows the TeV counterpart of the X-ray flare. The chance probability of the largest 
HAWC flux measurment to occur as a fluctuation above the constant average flux, considering a 
trial factor for 56 time bins, is 1.9 • 10 4 (3.7a). A second X-ray flaring state was observed in 
May 2014 in the lower 0.3-10 keV energy band by the X-ray Telescope (XRT) onboard the Swift 
satellite [p4] and is included in Fig. [] as the vertical blue band on the right side. No clear flux 


6 




















































Monitoring TeV blazars with HAWC 


Robert J. Lauer 


feature is observed in the HAWC data around that time. 

The green vertical band indicates a period where the VHE gamma-ray telescope FACT observed 
several flares from Mrk 501 with flux variations changing on the time scale of hours [21]. This 
triggered a follow-up observation by the VHE gamma-ray telescope H.E.S.S. [j25l], with a peak flux 
at MJD 56832. None of the weekly HAWC flux measurments during that period are larger than 2a 
above the flux average. It is important to note the transits time windows of the HAWC observations 
never overlap with the night time observations of either FACT or H.E.S.S. and that the averaging 
over 7 days in the HAWC analysis reduces the sensitivity to flux variability over less than one 
hour reported for this flare. The HAWC results presented here reveal that the week-averaged fluxes 
during these three weeks were not significantly higher than the average value over the whole year. 


4. Conclusion 

Data taken between June 2013 and July 2014 with an early configuration of the HAWC Ob¬ 
servatory, operating with approximately one third of the completed array, have been analyzed with 
a likelihood technique to produce flux measurments for 56 week-long time intervals. These light 
curves for the two clearly detected extra-galactic sources, BL Lacertae objects Markarian 421 and 
Markarian 501, both show indications of variability and both have the highest VHE flux values 
coinciding with X-ray flare observations. Detailed joint analyses of the multi-wavelength data and 
a finer time binning for selected periods of interest could reveal more features of the VHE gamma- 
ray flux during known flares of these sources. 

New HAWC data with improved sensitivity from several months of operating a grown configura¬ 
tion and, after March 2015, the completed HAWC Observatory are now available for analysis. With 
the increased sensitivity of the larger array, a light curve analysis with reduced flux uncertainties 
and shorter interval lengths of 1 transit will be possible and can provide regular monitoring for any 
blazar, or in fact any sky location, in the field of view of HAWC. 
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